Fidelity susceptibility in the two-dimensional transverse field Ising and XXZ models 
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We study the fidelity susceptibility in the two-dimensional(2D) transverse field Ising model and the 
2D XXZ model numerically. It is found that in both models, the fidelity susceptibility as a function 
of the driving parameter diverges at the critical points. The validity of the fidelity susceptibility to 
signal for the quantum phase transition is thus verified in these two models. We also compare the 
scaling behavior of the extremum of the fidelity susceptibility to that of the second derivative of the 
ground state energy. From those results, the theoretical argument that fidelity susceptibility is a 
more sensitive seeker for a second order quantum phase transition is also testified in the two models 
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I. INTRODUCTION 

Fidelity, a concept emerging from quantum informa- 
tion theory, has recently become an attractive approach 
towards the study of critical phenomena in condensed 
matter physics. In a quantum many-body system, the 
quantum phase transition is completely driven by the 
quantum fluctuation in the ground state and is incar- 
nated by an abrupt change in the qualitative structure 
of the ground state wavefunction as the system varies 
across the critical point 1]. Therefore, being a measure 
of the similarity between two states, the fidelity is ex- 
cepted to show a dramatic change across the transition 
points. This motivated people to start exploring its role 
played in quantum phase transitions [3, H, Q- More- 
over, as the fidelity can be viewed as a space geometrical 
quantity, no a priori knowledge of the order parameter 
and symmetry breaking of the system is required. This 
is thus a great advantage to the study of quantum phase 
transitions using fidelity approaches. 

Following the streamline of fidelity, some alternative 
schemes, like the fidelity susceptibility fidelity per 
site [^.operator fidelity and density- functional fi- 
delity 3, have been proposed. As to establish a closer 
picture to condensed matter physics, we follow the con- 
cept of fidelity susceptibility in this paper. Mathemati- 
cally, the fidelity susceptibility is just the leading term of 
the fidelity. It defines the response of the fidelity to the 
driving parameter. As a result, the singularity of the fi- 
delity across the transition points could thus be reflected 
in the divergence of the fidelity susceptibility. In fact, 
this argument has been consolidated by the results in a 
number of one-dimensional quantum many-body systems 
(See also a review article flol|). 

In this paper, we investigate the behavior of the fi- 
delity susceptibility in two two-dimensional (2D) models. 



namely the 2D transverse field Ising model and the XXZ 
model numerically. Our results show that the fidelity 
susceptibility as a function of the driving parameter di- 
verges at the quantum phase transition points in both of 
the models. The scaling behavior of the extremum of the 
fidelity susceptibility at the transition point and that of 
the second derivative of the ground state energy are also 
compared. From those results, the theoretical argument 
that the fidelity susceptibility is a more sensitive indica- 
tor than the second derivative of the ground state energy 
in searching for a second order quantum phase transition 
is testified. Besides, it is also found that the fidelity sus- 
ceptibility shows a scaling behavior in the vicinity of the 
critical point and its critical exponents for both models 
are also obtained through finite-size scaling analysis. 



II. FORMULISM 



For a general form of the Hamiltonian, 



(1) 



where Hj is the driving Hamiltonian and A denotes its 
strength. The fidelity is the modulus of the overlap be- 
tween two ground states |5'o(A)) and |\I'o(A + 5\)) Q, 



F(\,\ + 5\) = |(*o(A)|*o(A + <5A))|. 



(2) 



Since our focus is on continuous quantum phase tran- 
sitions, the ground state of the Hamiltonian is non- 
degenerated for a finite system. |^o(A + ^A)) can thus 
be obtained from the time-independent non-degenerated 
perturbation theory. Extracting the leading term of the 
fidelity, the fidelity susceptibility can be expressed as 



E 



|(^„(A)|g,|^o(A))|^ 
[£;„(A)-i?o(A)]2 
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where |\E'„(A)) is a set of orthogonal basis satisfying 
i7(A)|*„(A)) = i;„(A)|*„(A)). On the other hand, con- 
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FIG. 1: (Color online) Two-dimensional structures for difler- 
ent system sizes A*' = 10, 16, 18, 20 which can be placed on a 
square lattice with periodic boundary conditions. 



sider the second derivative of the ground state energy 
with respect to A jl2| . 
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(4) 



one may easily realize that the above expression is very 
similar to the perturbation form of fidelity susceptibihty 
in Eq. ^ except having different exponent in the denom- 
inator. Therefore, one may expect that both the singular- 
ity of fidelity susceptibility and the second derivative of 
the ground state energy arc intrinsically due to the van- 
ishing of the energy gap in the thermodynamic limit . 
However, the difference in the exponent of the dominator 
makes fidelity susceptibility a more sensitive quantity in 
searching for quantum phase transitions. That is to say, 
while the fidelity susceptibility shows a divergence at the 
critical point, the second derivative of the ground state 
energy may still be a continuous function. 

Furthermore, to study the scaling behavior of the fi- 
delity susceptibility around the critical point, we may 
perform finite-size scaling analysis [H, Let's con- 

sider a system consisting of N sites such that N ~ L"^, 
where d is the real dimension of the system. Around the 
critical point Ac, the fidelity susceptibility behaves as 



Xf(A) 



1 



|Ac-A|"^ 



(5) 



where a"*" (a ) is the critical exponent of the fidelity sus- 
ceptibility above (below) the critical point, is the 
quantum adiabatic dimension and hence xfW/L'^'' is 
an intensive quantity. For a finite system, if the fidelity 
susceptibility shows a peak at a certain point Amax, it's 



maximum value scales like 

XF(Amax) i'*", 



(6) 



where d'^ is the critical adiabatic dimension. The above 
two asymptotic behaviors satisfy [l3| 



Xf(A,£) 



A 



L-'ia+da +B(A- An 



(7) 



where A is a constant, B is a non-zero function of A and 
both of them are independent of the system size. From 
Eq. ([7]), one can find that the rescaled fidelity suscepti- 
bility is an universal function of L''(A — Amax), 



XF(Aniax,-^) - Xf(A,L) 

Xf(A,L) 



= /[L'^(A-A„,ax)], (8) 



where v is the critical exponent of the correlation length. 
The critical exponent of the fidelity susceptibility can 
then be obtained as [3, fisj 



a a 



(9) 



III. TWO-DIMENSIONAL TRANSVERSE 
FIELD ISING MODEL 

The Hamiltonian of the 2D transverse field Ising model 
[T^ [TtI defined on a square lattice reads 



-^Ising — 



h 



(ij) 



(10) 



where Sf, Sf , S?{S^ = /2, k = x,y,z) are spin-1/2 
operators at site i, h is the transverse field strength in 
unit of the Ising coupling, and the sum () runs over 
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FIG. 3: (Color online) The second derivative of the averaged 
ground state energy of the 2D transverse field Ising model 
on a square lattice as a function of h. The inset shows the 
scaling behavior of the minimum of the second derivative of 
the averaged ground state energy. The slope of the line is 
approximately 0.103. 



the nearest-neighboring pairs on the lattice. Periodic 
boundary conditions are assumed. This model was orig- 
inally introduced by de Gennes to describe potassium- 
dihydrogen-phosphate type ferroelectrics [ig] and has 
been studied extensively via various approaches, like real- 
space renormalization group [T^, density- matrix renor- 
malization numerical diagonalization[20l [2l|, and 

entanglement [22| . 

Obviously, the Hamiltonian of the model commutes 
with the parity operator = Hi ■ Thus each eigen- 
state of the Hamiltonian is also an eigenstate of P. The 
Hilbert space can then be decomposed into subspace V{p) 
where p is the eigenvalue of P and is specified in each 
subspace. For a finite system, the ground state of the 
Ising model is non-degenerate in each subspace, thus the 
perturbation expansion as introduced in the previous sec- 
tion is valid as long as the lattice is finite. In the ther- 
modynamic limit, the model exhibits a quantum phase 
transition at 1/hc ~ 0.328 [13, HI. For h > h^, the 
transverse field dominates and the ground state is a para- 
magnetic phase, with spins almost fully polarized in the 
z-direction. For h <gC he, the ground state is a ferromag- 
netic phase and is doubly degenerated. 

To study a model on a two-dimensional square lattice 
with periodic boundary conditions, we need to construct 
proper lattice structures that are suitable for exact diago- 
nalization. In this paper, we will diagonalize two models 
with system sizes N = 10, 16, 18, 20, whose structures 
are shown in Fig. [T] The effective length L — \/N might 
then be a real number instead of an integer. 

Fig. [2]shows the numerical result of the fidelity suscep- 
tibility of the Ising model on a square lattice for various 
system sizes. It can be seen that on both sides around 
the critical point, the averaged fidelity susceptibility is 
an intensive quantity, i.e. xf ^ N ^ L'^ and we have 
= 2. More importantly, the averaged fidelity suscep- 
tibility for different N all show a peak at /imax- This peak 



position of the fidelity susceptibility /imax is plotted as a 
function of 1/iV, as shown in the inset of Fig. O The 
linear fitting gives 



2.95 



6.56 



In the thermodynamic limit, we obtain 

K = 2.95 ±0.01, 

which gives 1/hc = 0.326 ± 0.001. Comparing this value 
to the critical point 1/hc — 0.328 obtained in previous 
works [13, HH , our result here is consistent with them up 
to two digits. 

Moreover, we can also see that the averaged fidelity 
susceptibility peaks sharper for a larger N and is in fact 
scales approximately with iV°-^^ as shown in the inset of 
Fig. [4] Physically, as the ground state wavefunction of 
the model changes abruptly across the transition point, 
the fidelity susceptibility, as a measure of the leading 
response of the fidelity to the driving parameter, is in- 
tuitively expected to show a divergence at the critical 
point. Here, we have shown numerically this is in fact 
the case and verified the significance of fidelity suscep- 
tibility in signaling for the quantum phase transition in 
the 2D Ising model. 

Fig. [3] shows the second derivative of the averaged 
ground state energy of the Ising model for several sys- 
tem size of a square lattice as a function ol h. As it is 
well-known that the Ising model exhibits a second order 
phase transition, the second derivative of the averaged 
ground state energy is expected to show a minimum at 
the transition point. From the inset of Fig. [3l it is found 
that the minimum value of the second derivative of the 
averaged ground state energy scales approximately with 
^0.103^ Comparing this value with that of the fidelity 
susceptibility, which is about 0.51 (d^ ~ 3.02), we may 
conclude that the fidelity susceptibility is a more sensi- 
tive tool in detecting for a second order quantum phase 
transition. 

Fig. [4] shows the finite-size scaling analysis in the case 
of power-law divergence of the 2D transverse field Ising 
model. The rescaled fidelity susceptibility collapsed to 
a single curve for various system sizes. The critical ex- 
ponent of the correlation length can thus be obtained as 
1/ ~ 1.40. Together with the slope of the line in the inset 
of Fig. [Hand from Eq. ([5]), the critical exponent of the 
fidelity susceptibility is found to be 



1.02 

L40 



0.73. 



IV. TWO-DIMENSIONAL XXZ MODEL 



For the 2D XXZ model, the Hamiltonian is given by 



Hxxz 



(ij> 



{s?s? + sysy + As?s?), 



(11) 
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FIG. 4; (Color online) The finite-size scaling analysis is per- 
formed for the case of power-law divergence for the 2D trans- 
verse field Ising model. The fidelity susceptibility for different 
system sizes is a function of N"^'^{h — /imax) only, with the 
critical exponent v ~ 1.40. The insert shows the scaling be- 
havior of the maximum of the averaged fidelity susceptibility. 
The slope of the line is 0.51. 
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FIG. 5: (Color online) The averaged fidelity susceptibility in 
the ground state of the 2D XXZ model on a square lattice as 
a function of A. The inset shows Amax as a function of 1/N . 
The i/-intercept of the line is 1.05 ± 0.02. 
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FIG. 6: (Color online) The second derivative of the averaged 
ground state energy of the 2D XXZ model as a function of A. 
The inset shows the scaling behavior of the minimum of the 
second derivative of the averaged ground state energy. The 
slope of the line is approximately 0.96. 




FIG. 7: (Color online) The finite-size scaling analysis is per- 
formed for the case of power-law divergence for the 2D XXZ 
model. The fidelity susceptibility for iV > 10 is a function of 
N"^'^(^^ — Amax) only, with the critical exponent v ~ 3.00. 
The insert shows the scaling behavior of the maximum of the 
averaged fidelity susceptibility. The slope of the line is 1.81. 



where A = Jz/ Jy{Jx = Jy) is the dimensionless parame- 
ter characterizing the anisotropy of the model. The sum 
is over all nearest-neighbors on a square lattice. Again, 
periodic boundary conditions are assumed. For the XXZ 
model in two dimensions, there exists no exact solution. 
One has to use either approximate analytical approach 
such as the spin- wave theory or numerical approach such 
as exact diagonalization studies of a finite lattice. For 
the latter approach, to obtain results in the thermo- 
dynamic limit, finite-size scaling analysis must be per- 
formed 26, 27]. Therefore, a physical quantity which is 
more sensitive to the system size than the traditional sec- 
ond derivative of the ground state energy would be very 
useful to study the critical phenomena numerically. 

From Eq. (fTTj) , it can be easily seen that the Hamilto- 
nian of the XXZ model commutes with the z-component 



of total spin operator ^^otai = Si ^i- Thus, each eigen- 
state of the Hamiltonian is also an egienstate of 
The Hilbert space can then be decomposed into numer- 
ous subspaces V{M), where M is the eigenvalue of -S'^otai- 
For a finite sample, the ground state of the XXZ model is 
non-degenerated in any of the admissible subspace V{M) 
[i^li^l. Therefore, the perturbation expansion can also 
be applied to this model as long as the system is finite. In 
the thermodynamic limit, the quantum phase transition 
takes place at the isotropic point Ac — 1. This phe- 
nomenon can be understood by the picture of the first 
excited energy levels crossing at the transition point [25j . 
For A ^ Ac, the last term in the Hamiltonian dominates 
and the ground state is an antiferromagnetic phase along 
the z-direction. For A <C Ac, the first two terms in the 
Hamiltonian dominate and the ground state is also an 
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antiferromagnetic phase, but in the xy plane. It is well 
known that long-range orders are present in both of the 
two phases. However, whether there's a long-range order 
at the critical point is still an open question. With the 
help of fidelity susceptibility, we may also find some hints 
towards this question. 

The numerical result of the averaged fidelity suscepti- 
bility for various system sizes of the 2D XXZ model on 
a square lattice as a function of A is shown in Fig. [5] 
The averaged fidelity susceptibility is an intensive quan- 
tity, meaning that xf ^ N, on both sides of the criti- 
cal point. Moreover, like the previous case of the Ising 
model, the averaged fidelity susceptibility of the XXZ 
model also shows a peak at Amax- The inset of Fig. [5] 
shows the peak position of the fidelity susceptibility Amax 
as a function of 1/N. The linear fitting gives 



In the thermodynamic limit, we obtain 

Ac = 1.05 ± 0.02. 

Comparing to the theoretical critical point Ac = 1, our 
result here is consistent up to two digits. 

Besides, from the slope of the straight line in the inset 
of Fig. [3 it is found that the peak of the averaged fi- 
delity susceptibility scales with the system size like N^-^^. 
Therefore, one may expect the fidelity susceptibility to 
show a singularity at the critical point in the thermody- 
namic limit. Hence, the validity of the fidelity suscep- 
tibility as a seeker for the quantum phase transition is 
also verified in the 2D XXZ model. Nevertheless, follow- 
ing the idea of the implication of existence of long-range 
correlation from the divergence of the fidelity suscepti- 
bility [2^ , we argue that long-range correlation is in fact 
present at the transition point of the 2D XXZ model. 
This is also in agreement with the previous conclusion 
drawn from the study of the 2D XXZ model using entan- 
glement [29]. 

In comparison, the second derivative of the averaged 
ground state energy for various system sizes exhibits a 



minimum at the transition point, as shown in Fig. [6l 
From the inset of Fig. [51 it is also found that the min- 
imum value of the second derivative of the averaged 
ground state energy scales approximately with iV°-^^, 
meaning that it shows a slower divergence at the criti- 
cal point compared to the fidelity susceptibility. In other 
words, the fidelity susceptibility is again a more sensitive 
candidate in seeking for the quantum phase transition in 
the 2D XXZ model. 

Fig. [7] shows the finite-size scaling analysis in the case 
of power-law divergence of the 2D XXZ model. The 
rescaled fidelity susceptibility almost collapsed to a sin- 
gle curve for a large enough system size, say N > 10. 
The exponent of the correlation lenght is obtained as 
1/ ~ 3.00. From the slope of the inset in Fig. [7] and using 
Eq. © , the critical exponent of the fidelity susceptibility 
is calculated to be 

3.62 

a = = 1.21. 

3.00 



V. SUMMARY 

To conclude, through the numerical study of the fi- 
delity susceptibility in the 2D transverse field Ising model 
and the 2D XXZ model, we found that the fidelity sus- 
ceptibility as a function of the driving parameter diverges 
in both models at the critical point. By comparing the 
scaling behavior of the extremum of the fidelity suscep- 
tibility to that of the second derivative of the ground 
state energy, we also showed that fidelity susceptibility 
is a more sensitive indicator in detecting for a second or- 
der quantum phase transition. By performing finite-size 
scaling analysis, the critical exponent of the fidelity sus- 
ceptibility in both models are also obtained. Finally, due 
to the divergence of fidelity susceptibility in the 2D XXZ 
model, we argued that the system shows a long-range 
correlation at the critical point of the model. 
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China (Project No. CUHK 400807). 
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